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Abstract We describe a novel analogue algorithm that allows the simultaneous fac¬ 
torization of an exponential number of large integers with a polynomial number of 
experimental runs. It is the interference-induced periodicity of “factoring” interfero- 
grams measured at the output of an analogue computer that allows the selection of 
the factors of each integer mi2ii3iHi . At the present stage the algorithm manifests an 
exponential scaling which may be overcome by an extension of this method to corre¬ 
lated qubits emerging from n-order quantum correlations measurements. We describe 
the conditions for a generic physical system to compute such an analogue algorithm. 
A particular example given by an “optical computer” based on optical interference 
will be addressed in the second paper of this series Q. 

Keywords quantum computation • interference • algorithms • analogue computers • 
factorization • exponential sums • Gauss sums 


1 Fundamental principle of our work 

Multiplying numbers is much easier than the inverse problem of finding the factors 
of large integers. Indeed, the security of codes relies on the current inability of a 
fast solution of this problem but is endangered by the well-known Shor’s factoring 
algorithm Il6ll7l[^l9lfra . 

Recently, important works about the use of different kinds of exponential sums 
ifTTl for factorization purposes have been published 012lll3lfT4lfT5lfT6lfT7lfT8lll9ll20l 
l2ni22l . Our method differs from the past experimental realizations Il^l24ll25Il26ll27l 
I28ll29l in three important simultaneous achievements |[il|2l[3 . First, the division of N 
by the test factors i is not pre-calculated, but it is performed by the experiment itself. 
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Second, several test factors are tested simultaneously. Third, a scaling property inher¬ 
ent in the recorded interferograms allows us to obtain the factors of an exponential 
number of large integers. 

The core of Shor’s factoring algorithm stands on the measurement of the period¬ 
icity in the dominant maxima of the quantum probability distribution at the output 
of a quantum computer ifTOll . In line with Shor’s idea, the key behind the algorithm 
we present in this paper stands in the measurement of the periodicity in the max¬ 
ima of Continuous Truncated Exponential Sums (CTES) by performing first-order 
“factoring” interference processes with a physical system. Interestingly, the number 
n of necessary experimental runs scales logarithmically with respect to the largest 
integer to be factored. A noteworthy theoretical result at the core of the algorithm 
is that the periodicity of the resulting n interference patterns as a function of a con¬ 
tinuous physical parameter in a given range, when appropriately scaled, allows us to 
achieve factorization of large numbers by simply looking at the interference maxima 
at integer values. 

The paper is organized in the following way. In section II will be given some 
mathematical background about CTES in connection with the hyperbolic function. 
In section III we will describe the factoring algorithm for a generic physical system 
and the conditions this system must satisfy. We will provide a generalization of the 
described algorithm in Section IV. Section V and VI will address final remarks and 
perspectives of extensions to polynomial scaling methods of factorization, respec¬ 
tively. 


2 Hyperbolic function, CTES periodicity and factorization 


We define a continuous truncated exponential sum (CTES) in the form |[il|2l[3|4l 

= (1) 

where is the modulo-squared value of the generalized curlicue function lISOl fTl 

1 




1 ^ 

— £ exp[27r/(m-l)^C] 


2 


( 2 ) 


of integer order j> I, with M >2 interfering terms, and 

( 3 ) 

is the hyperbolic function with continuous variable 0 <^< 1 . 

In Eig. we represent the modulo squared of the curlicue function = 

5 (^ 1 !) in its dependence on the argument for M = 3,4,5 and y = 1,2,3. We note 
that p has, for all the orders j, a dominant maximum at 1 ^ = 0 with ( 0 ) = 1 
and decaying oscillations on the sides. The higher the order j and the truncation pa¬ 
rameter M of the curlicue function are, the larger is the number of decaying oscil¬ 
lations on the sides and the sharper is the dominant peak. Moreover, we recognize 
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from Eq. (|^ the periodicity property (i^ -(- 1) = Therefore, it suffices 

to consider in the domain—1/2< ^ < 1/2. In addition, {Q\^ 

is symmetric with respect to = 0. 



Fig. 1 Modulo-squared value of the generalized curlicue function in its dependence on the 
argument ^ for the number M = 3,4,5 of interfering waves (column) and the power j — 1,2,3 of 
the phase shift (row). For increasing M the dominant peak becomes nan'ower, which will make 
it easier in our algorithm to check if a peak corresponds to a factor or not. Unfortunately, at the 
same time the number of side maxima increases as well. However, we note that for a fixed j and 
increasing M the value of the maxima of second order decreases and becomes sharper. On 
the other hand, for a fixed M but increasing j the value of the maxima of second order increases 
and becomes wider (4j- 


The hyperbolic function / in Eq. ^ induces in the function a notable 

periodicity. Indeed, the function is characterized by dominant maxima, which 

repeat each time / assumes integer values. 

Why does such a CTES periodicity matter in factorization? 

In order to answer this question we first point out that, as shown in Refs. |[T] 
Em, the factorization problem could be, in principle, solved if we can achieve the 
complete knowledge of the hyperbolic function /. 

Indeed, if we look at / as a function of the new variable obtained by the scaling 
relation 




( 4 ) 
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we obtain 



(5) 


For each possible value of N, the factors are given by the integer values = £ such 
that 


N 

m = j=k, 


(6) 


with k a positive integer. 

Unfortunately, it is not an easy task to compute the hyperbolic function so that for 
any given integer N the condition (j^ can be verified in order to identify the factors. 

Interestingly, we can exploit the constructive/destructive periodic interference 
characterizing the CTES function in Eq. 0 as a tool in the distinction between fac¬ 
tors and non factors of any given number N. In particular, the rescaled hyperbolic 
function f{^N) in Eq. 0 corresponds to the rescaled CTES 


2 


1 “ 

= — £ exp [2ni{m- l)V(‘?iv)] 

m=\ 


(7) 


as a function of = . The condition 0 leads to total constructive interference in 

the rescaled CTES function in Eq. 0. Indeed, the factors of an arbitrary 

number N are the integer values £ of corresponding to dominant maxima in the 
function 

In Eig.j^is simulated the rescaled CTES function in Eq. 0, with M = 3, j = 2, as 
a function of the variable G [330.84,337.21] for the factorization of = 111547. 
We can see that the two factors £ = 331,337 (represented by stars) correspond to 
complete constructive interference. On the other hand, for the other test factors (rep¬ 
resented by triangles) there is partially destructive interference. Moreover, there are 
absolute maxima (represented by points) which do not correspond to integer test fac¬ 
tors. 

In Eig. instead, we have simulated the CTES function in the case of j = 3 for 
the same value of N and M and the same range of values of ^N- As expected, it turns 
out that as the order j of the exponential sum increases, the peaks associated with the 
absolute maxima become sharper. On the other hand, the values of the second order 
maxima in the interference pattern increase for larger orders j. In order to suppress 
such maxima, it is necessary to increase the number of terms M in the sum. Also, in 
such a case, the first order peaks become increasingly sharper. 

If we consider the range of values of in Eig.j^and Eig.[^ the closer a non factor 
is to a factor, the larger the corresponding value is of the rescaled CTES function. Of 
course, if we move to integer values of farther from the factors there is a larger 
probability of partial or total destructive interference between the terms in the CTES 
function. This makes the distinction between factors and non factors easier. 

In general, the maximum possible step A ^ between consecutive values of ^ where 
the CTES has to be computed can not be exactly determined since it depends on the 
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Fig. 2 Rescaled CTES (j=2) function (^jv) in Eq. with M = 3 for = 111547, as a 

function of the variable ^ in the interval [330.84,337.21] (3. We can see that the two fac¬ 
tors — 331,337, represented by stars, correspond to complete constnictive interference, with 
respect to the other integer trial factors — 332,333,334,335,336, represented by triangles, 
which present partially destructive interference. 


value of the factors of the integer N < NAAmax to be factored. However, it can be 
easily shown m that such value is included in the interval 

N-^<A^<1 ( 8 ) 


if G or the interval 


V]^<A^<1 (9) 

if^NG[VN,N]. 

In conclusion, the CTES function in Eq. o allows us to extract the information 
about factors encoded by its periodicity in the dominant maxima. Such a periodicity 
is imprinted by its functional dependence on the hyperbolic function f{^) = I 
We recognize all the values ^ corresponding to an integer value of /(^) = l/^ as 
dominant maxima in the interference pattern. When for one of these values of ^ we 
find, for a given large number N, be an integer, such an integer is a factor 

ofV. 

However, a digital implementation of the CTES function for factoring purposes 
would rely on performing an exponential number of divisions, which are expensive 
operations for digital computers. On the other hand, an analogue implementation 
of such a “factoring” function would rely on the help of “nature” to perform such 
divisions for us. In the next section, we will describe the prerequisites for a generic 
physical system to compute a factoring algorithm based on the implementation of 
CTES interferograms. 
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Fig. 3 Rescaled CTES (j=3) function in Eq. Q with M — 3 for V — 111547, as 

a function of the variable ^ in the interval [330.84,337.21] (2). The two factors — 
331,337, represented by stars, correspond to complete constructive interference, with respect to 
the other integer trial factors — 332,333,334,335,336 in such an interval, represented by 
triangles. As expected, the peaks associated with the absolute maxima in the case 7 = 3 are 
sharper than the respective peaks in the case j — 2 represented in Fig.[^ On the other hand, the 
value of the maxima of second order in the function increases at the increasing of the 

order j. 


3 Factoring analogue algorithm 

In the previous section we have shown how the implementation of the CTES func¬ 
tion in Eq. 0 would allow, in principle, the prime number decomposi¬ 

tion of arbitrary integers. Now, we want to describe the analogue implementation of 
the CTES algorithm with a generic physical system that exploits interference. This 
will lead us to the introduction of a two-dimensional interferogram = 

as a function of a continuous physical parameter and a discrete 

physical parameter x associated with two independent observables and Ox, respec¬ 
tively. In the second paper of this series 0 we will give an example of an “optical 
computer” based on a multi-path Mach Zehnder interferometer, where the physical 
parameter and x will correspond to the wavelengths associated with the spectrum 
of the optical source and the optical-path unit, respectively. 

We will first introduce a CTES “factoring” procedure which takes advantage of 
a single interferogram (o^ ;x) associated with a given value of x. We will deter¬ 
mine the range Nmin < N < N^ax of factorable numbers N by covering all the trial 
factors in either the range [1, yV] or [VV, V] with a given range o,nin 1^o^< Omax of 
values of the observable m. 

Next, we will show that the CTES “factoring” algorithm can exploit the degree of 
freedom in the variable x of the analogue function in order to check all 

the necessary trial factors of an exponential number of large integers N m. In par- 
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Fig. 4 Rescaled CTES {j = 2) function in Eq. with M = 3 for M = 111547, 

as a function of the variable = in the interval [230.9,237.1], We can clearly see that 
all the integer trial factors = 231,232,233,234,235,236,237 in such a range, represented 
by triangles, correspond to a relatively small value of the CTES function with respect to the 
dominant maxima so they can be easily disregarded as possible factors. 


ticular, such a procedure is based on the measurement of the periodicity of a number 
n of interferograms for different suitable values of x. We will interestingly hnd that 
the value of n scales logarithmically with respect to the largest number Nmca to be 
factored. 


3.1 Analogue implementation of the CTES function 

We have shown that the implementation of a CTES function in Eq. Q would allow 
the factorization, in principle, of arbitrary numbers. Unfortunately, the calculation 
of such a sum would require an exponential number of divisions associated with the 
computation of the function /(^). On a digital computer, for which division is a rather 
costly process, such a computation turns out to be very slow. Therefore, it would be 
more efficient to reproduce the CTES function with an analogue technique in order to 
solve the problem quickly. More explicitly, an analogue implementation of the CTES 
algorithm is possible if there is a physical system able to compute divisions for us 
and then read out the factors by taking advantage of the interference process leading 
to CTES interferograms. 

In particular such a physical system needs to fulhll three main requirements. 

Eirst, we assume that the system is characterized by two independent observables 
Ox and . The observable Ox can be tuned to values 


( 10 ) 
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with j a positive integer and m= 1, The unit of measurement x can be, in prin¬ 
ciple, arbitrarily varied. On the other hand, the observable assumes a continuous 
range of values 


o^=^x, ( 11 ) 

where x is the unit of measurement chosen in Eq. ( [T0| ). 

Second, the physical system needs to be able to compute the hyperbolic function 
in Eq. (|^ in the form of ratios between the values of Ox and for a given interval 
of the variable ^. 

Third, the system must be able to exploit interference in order to reproduce the 
interferogram 




1 M 


m=l 


.Ox 


mj) 






( 12 ) 


as a function of = ^x. Of course, for a parallel evaluation of the sum for several 
values of ^ it is necessary that the system contains at the same time the information 
about all the possible values of the physical observable . The corresponding CTES 
function in Eq. 0 will finally allow us to extract the information about 

factors. 

In particular, for a generic value of the number N to be factored, it is possible to 


look at the obtained interferogram in Eq. (12 1 as a function of the rescaled variable 
in Eq. 0 


=N^ = —O^, 


(13) 


where we use Eq. (11 1 in the second equality. 


Indeed, we obtain the rescaled interferogram 




(14) 


corresponding to the rescaled CTES function defined in Eq. (|^. 

Each time there is a dominant maximum at a value of for which in Eq. (131 
is an integer, we find such an integer to be a factor of N. 

It is important to point out that the information about the one-dimensional CTES 
function ) in Eq. 0 is inferred by the two-dimensional CTES interferogram 

in Eq. Jl2|l as a function of the physical variable and the physical 
parameter x. Indeed, asnecomes clear later, this feature will turn out to be one of the 
key points to understand the working principle behind the CTES analogue algorithm. 
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3.2 Factorization with a single interferogram 


In this section, we describe the “factoring” procedure based on the use of a single 
interferogram given by the CTES analogue function /(^’^) {o^;x) in Eq. Jl2b recorded 
at a given value of x. We address the question of the interval Nmin,x < Nmax x of 
numbers N factorable by covering all the trial factors in either the range [3, V^] □or 
[v/iV, A^] with a given range o„„„ <o^< o^ax of values for the observable . 

In general, for each integer N a generic trial factor i can be checked only if 

{: girnax)-, _ rN N , 

SA — € [sa/' 5 \ — I. ^ ^ ? (15) 


where and are respectively the smallest and largest values that the vari¬ 


able in Eq- < 13 I can assume for the rescaled interferogram {^}^\x). 

We consider the case in which we want to check all the trial factors £ S [3, v/A^] 


leading from Eq. (15 1 to the condition 


N 


N 


Method (1); = £ e [3, ViV] C 


(16) 


By dividing the upper and lower bounds of each interval respectively by ■\/N and 3 
we find 


which implies 


1 G 


N 



^ <-1 

Ojnin ^ Omax 

5x X 


(17) 


(18) 


for each integer N in the interval Nmin,x < Nmaxjc to be determined. By squaring 
the third term in this series of inequalities we obtain the condition 


x< X 


( 1 ) = 


3o^ 

^min 


(19) 


giving an upper bound to the choice of the parameter x associated with the rescaled 
interferogram /(^’D {^;x) independent of the number N to factorize. Erom Eq. (181 
we also easily obtainR 


= 

^^max,x — 


3x 


and 


= 

min.x 


r 1 


( 20 ) 


* We are sure that there is at least one factor of N in such interval. The trial factor 2 is obviously 
excluded since it is easy to recognize if N is an even integer. 

^ We recall that for any real number y the ceiling \y] is the smallest integer lai'ger than jc, while the floor 
[yj is the largest integer lower than y. 
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defining the interval of factorable numbers with a single interferogram associated 
with a given value x in Eq. (19 1 in the range o,nm ^ < Omax- 

In the particular case of an interferogram recorded at the maximum possible value 
X = in Eq. (19 1 , the condition (18 1 reads 

N VN 

If Omax/omin is an integer, we find the largest but also the only integer 

q„2 

iV(i) = (21) 


factorable by using a single experimental interferogram {o^;x) in Eq. (12 1 recorded 
for the value x = x*. 

We consider now the second method in which we want to check all the trial factors 
= (■ & leading from Eq. to the condition 


/— N N 

Method (2); ^ ^ [v ,A^] ^ -^Omax]' 


( 22 ) 


By dividing the lower bound and the upper bound of each interval by v/iV and N, 
respectively, we find 


, _ r Ofnax n 

1 e - 0„ 


X X 


(23) 


which implies 


, ^ '\/N ^max ^ , 

1 > - 0,nin < - > 1 

X X 


(24) 


for each integer N in the interval Nmin,x <N< Nmaxx to be determined. Erom the last 
inequality in Eq. ( p4j l we obtain the condition 

X<X^^'^ =Omax- (25) 


Erom Eq. (24 1 we also easily obtain 


= 

^^max,x — 


L ^min J 


and 
= 1 

-'’mm — 


(26) 


defining the interval of factorable numbers with a single interferogram associated 
with the generic value x in Eq. (25 1 for the given range o„„>, < < Omax- 

We consider now the case oFa single interferogram recorded at the maximum 
value X = x^^') = Omax in Eq- <25 i leading to the largest interval 


Wp) = 1 < W < = 

min — ^ _ -I'' _ ^^max — 


max 

^2 


( 27 ) 
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of factorable integers. 

We finally demonstrate that with a single interferogram it is possible to factorize 
a number 

^min 


of integers exponential with respect to the number of binary digits rimax associeted 
with Nifuix- However, in general, the largest factorable integer Nmax is limited by the 
value Omax/omin associated with the physical system. For this reason, in the next sec¬ 
tion we will describe a factorization procedure which takes advantage of several in- 
terferograms (o^ ;x) in Eq. (121 at different values x in order to factor numbers 
exploiting a limited fixed range oTvalues of a physical observable . In such a 
method, the maximum factorable number Nmax will depend only on the largest value 
achievable for the parameter x. 


3.3 Factorization with a sequence of interferograms 


So far we have restricted ourselves to a factorization method involving a single inter¬ 
ferogram /(^’A (o^;x) in Eq. (|l2| defined at a fixed value of the parameter .r. However, 


the remarkable scaling property = No^ jx characterizing the function (o^ ;x) 
allows us to vary the physical values of both and x. This implies that we can, in 
principle, change arbitrarily both the number we are looking at and the relative trial 
factors to check by simply varying these two physical “knobs”. 

In particular, we seek to address the question of the largest number Nmax that 
can be factored by measuring several interferograms I^^’-’\o^-,x) corresponding to 
different values of x in a hxed domain Omin Omax of the variable . 


3.3.1 Factorization of a single integer N 


We consider first the case of a single generic large integer N to be factored. The 
scaling property = No^ jx, written as At = ^xixjo^ , implies that, in principle, there 
is no limit to the largest possible value of N if we can vary the unit x up to arbitrary 
large values. Of course, in practice, the maximum achievable value of x is limited by 
the particular physical system we are using. We now demonstrate that the number n 
of values of x, for which the “factoring” interferogram fNl'l (o^; v) in Eq. (12 1 must 
be recorded, scales logarithmically with respect to the value of N. 

It should first be pointed out that a single interferogram registered at a given value 
X allows us to check only the trial factors 




^ (max) 




N 


— Cmin 1 Ofy 


(28) 


which, in general, for the fixed domain o„„„ < o^ < Omax of values o ^, may corre¬ 
spond only to a subset of the total range 3 <£< \/N or \/N <(.<N of values ^n = £ 
we would need to cover in order to factor a generic integer N. For this reason, we 






12 


Vincenzo Tamma 


consider a suitable sequence of values x = Xi, with i = 0,1, ...,n — 1. Each interfero- 
gram registered at the value x = x,, with / = 0,1, ...,n — 1, allows us from Eq. ( [28] l to 
cover all the trial factors 

N N 

^ ^ .i : .i+l] — { Otnim (29) 

Xi Xi 

with i = 0, l,...,n— 1, where 

N N 

^N.i+\ = Offuix = Oyfiifi (30) 

Xi Xi+\ 


satisfies the condition for consecutive intervals. 

This implies that the sequence x;, with i = 0, 1, associated with the n inter- 

ferograms is defined by the condition 


with 


and, from Eq. (29 1 , 


Erom Eq. (32 1 we obtain 


X,+ i = - < Xi, 

c 


,1-1-1 ^max ^ j 

^N.i Ofnin 


XQ ^ N 

Omin 


^N.n — 


leading to the number 


_ , ^N,n 

n = logcT — 
SV.O 


(31) 

(32) 

(33) 

(34) 

(35) 


of interferograms necessary to cover all the trial factors 

i: 0 ^ it. t 1 — i^^min ^Omin 

Xo Xo 

We deduce that the value xo for the first interferogram as well as the total number 
n of interferograms {o^;x) depend on the integer N to be factored and on the 
associated interval [^Arp ,^N,n] of trial factors to be checked. In particular, factorization 
can be achieved for a generic integer N only if 




Method (1); 

ie[3,VN]cg^^oAN,n] = [— 


X Ofyiifi 

xo 


(37) 
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or 


Method (2): 

= [\/N,N] C = [ 


(38) 


^ Ofnin ^ O n 


^0 


Xo 


c"]. 


Let us consider first the method (1). The lowest trial factor 3 to be checked leads 
to the condition 


^A ',0 = 


No„ 


Xo 


<3, 


(39) 


implying 


■^0 ^ 


N 


(40) 


In an analogous way, the value of the largest trial factor to be checked in Eq. (37 1 
determines the condition 


— 


No„ 


Xo 


c" > y/N. 


(41) 


This, together with Eq. (|40li, implies the minimum number 


= 

- 


< 


logc 


log. 


xo 


^min y/N 


y/N' 


N.min 


(42) 


of necessary interferograms (o^;x,), with i = 0,1,1, to factor a generic 
integer N with the method (1). 

We consider now the method (2) in Eq. (38 1 . The lowest trial factor y/N to be 
checked leads to the condition 


f- NOrnin , nr, 

^Nfi = - < VN, 

Xo 


which implies 


.( 2 ) 


^min Offiin 


(43) 


(44) 


In an analogous way, the value N of the largest trial factor to be checked in Eq. (381 
determines the condition 


f- ^^min n ^ at 

^N,n = - C" > N. 

XO 


( 45 ) 


This, together with Eq. (44 1 , leads to the minimum number 
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of necessary interferograms with i = 0,1, ...,n — 1, to factor a generic 

integer N with the method (2). 

This demonstrates that, in principle, the number n of interferograms necessary for 
factorizing an arbitrary large number N, using a given range [omin , Omm] of the physi¬ 
cal variable , scales logarithmically with respect to -s/N and thereby polynomially 
with respect to the number of binary digits associated with N. It is important to point 
out that the narrower the range Omin <o^< o^ax is, the closer the value of c in Eq. 


is to 1 and thereby the larger is the value of n in Eq. (35 i. However, it would be 
enough to cover a spectrum such that o^ax = '2-Omin in order to obtain a scaling that 
is logarithmic (base 2) and thereby polynomial with respect to the number of binary 
digits associated with N. 


3.3.2 Extension to the factorization of an exponential number of integers 


So far we have considered the case of factoring a single number with a sequence of 
n “factoring” interferograms. However, the scaling property = No^ jx implies, as 
pointed out before, that we can exploit the same interferograms in order to factor 
not only a single integer but any integer in any given range Nmm < N < Nmax- We 
determine how, in such a case, the number n of experimental tTins depends on the 
smallest and the largest number Nmin N^ax that can be factored. 

In particular, factorization can be achieved for all values of N, with Nmin <N< 
Nmax, only if for each single value is satisfied either the condition (37 1 for the method 


(1) or the condition (38i for the method (2). Let us consider first the method (1) in 
Eq. (0. The condition in Eq. (|3^ needs to be satisfied for each Nmin <N < Nmax, 
which implies 


NmaxOt 


max'xmin 


<3 


■*0 


leading to 


-^0 ^ -^0 _ Nmax 


(47) 


In an analogous way, the condition in Eq. (41 1 needs to hold for each Nmin <N< 
Nmax, leading to 


c"> 


1 


■^0 s/Nmin 

This, together with Eq. (|4^, implies the minimum number 


nxQ — 


> 


logc 


log, 


XQ 


n x/Nm 


OmmV ^'‘min 

N 


'ix/N~m 


mm 


(48) 


of necessary interferograms I^^'-l\o^;xi), with i = 0,1, to factor all the in¬ 
tegers N in any given interval Nmi,, <N< Nmax with the method (1). 
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We consider now the method (2) in Eq. (381. The condition (43 i needs to be 
satisfied for each Nmin < Nmax, which implies 


Omin ^ 1 

■*0 “ 


and thereby 


XQ 


.( 2 ) 


> 

Omin Ofnin 


— \/^m 


In an equivalent way, Eq. (50 1 reads 


nxQ — 


logc 


Xq 


> \togc\/Nmax\ — «l,(n 


(49) 


(50) 


of necessary interferograms I^^'-’\o^;xi), with i = 0,1, ...,n — 1, to factor all the in¬ 
tegers N in any given interval Nmm < Nmax with the method (2). 

In conclusion, the described algorithm allows in both method (1) and (2) the 
factorization of an exponential number 


AN — Nnuix Nm 


'Nm 


, '2J^max 


of integers, with timax number of binary digits of Nmax, by using a polynomial num- 

. The largest factorable 


mm ? '-'max 


ber of interferograms in a given physical domain [oj 
number is upper limited by the condition for xi^jornm in Efi- <47 1 and Eq. (49 1 . 


33.3 Example of factorization ofN < Nmax = 64 

We now describe the implementation of our factoring algorithm for Nmax = 64. We 
consider a generic observable with values o^ = ^x ina given range < o^ < 
Omax which satisfies the condition 


^ — Omax! Omin — 2 . 


(51) 


We consider first the method (1) in Eq. (371. It turns out that any integer in the 
range Nmin = 8 < < Nmax = 64 can be factored, according to Eq. (|48ll, by exploiting 


-„(i) 


log2 


64 

3V8 


= 3 


interferograms defined by Eq. (j^, with i = 1,2,3, 

where the value x = x^^^ associated with the first interferogram satisfies the condition 
in Eq. ( |47] l 

/V_ 64 

( 52 ) 


64 

T’ 
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From Eq. ([T3]l we can now determine all the values 


( 1 ) _ -i ( 1 ) 64 . 

X = X] ' r=C x}^' = —-I Omin, 


(53) 


with i = 0,1,2, for each of the n = 3 interferograms. We assume that our physical 
system is able to record such interferograms in the range Omin <o^< 2omin- In Fig- 
5 we simulate these interferograms, which are able to cover all the trial factors in the 
interval [3, v^] for any integer in the range A^„„„ = S <N < Nmax = 64. An example 
is given for the factorization of A = 15,63 by rescaling the axis as a function of 
^15 and ^63 according to Eq. (13i. The trial factors are marked, respectively, with 


continuous lines and dashed lines. We find the factors 3 and 5 of 15 corresponding to 
maxima of the interferogram On the other hand the factor 3 of 63 is 

associated with the maximum of the interferogram while the factors 

7 and 9 emerge from the maxima in the interferogram /(^’■')(o^ 

We consider now the method (2) in Eq. ( [3^ . According to Eq. (50 1 , any integer in 
the range Nmin = 1 < < Nmax = 64 can be factored by exploiting the same number 


^ = ^min- =3 


of interferograms I{o^-,xf''^) = (o^ ) defined by Eq. Jl^, with i = 1,2,3, 

.( 2 ) 


where the value x = Xq associated with the first interferogram satisfies the condition 


in Eq. (49i 


,,(2) 

^ = v^ = 8. 
Omin 


From Eq. 113 i we can now determine all the values 


A — — L Aq — o ' z, Umin i 


(54) 


(55) 


with i = 0,1,2, for each of the n = 3 interferograms. We assume that our physical 
system covers the spectrum Omin 2 o„„>, of values of a given observable . In 

Eig. 6 we simulate such interferograms, which are able to cover all the trial factors 
in the interval [a/A,^] for any integer in the range Am,„ = 1 < A < Nmax = 64. We 
again give an example for the factorization of A = 15,63 by rescaling the axis 
as a function of ^15 and ^^3 according to Eq. ( [q] ). We find that the factor 3 of 15 
corresponds to a maximum of the interferogram while the factor 5 is 

associated with a maximum of the interferogram On the other hand 

the factor 9 of 63 is associated with the maximum of the interferogram ). 
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Fig. 5 Factorization of any integer in the range Nmin = S <N < N,„a.x = 64 by using the method 
(1) aimed at checking all the trial factors in the rang e [3, VN]. We exploit n — 3 interferograms 
defined by Eq. jl^ , with i — 1,2,3, where the range of values 

satisfies the condition p 1 1 and the values are given by Eq. . We give an example for 
the factorization of A = 1 d 763 by rescaling the axis as a function or ^15 and <^63 according to 
Eq. (ig. The trial factors correspond, respectively, to continuous lines and dashed lines. We find 
the factors 3 and 5 of 15 corresponding to maxima of the interferogram On the 

other hand, the factor 3 of 63 is associated with the maxima of the interferogram 
while the factors 7 and 9 emerge from the maxima in the interferogram 


4 Generalization of the CTES algorithm 

In this section, we will introduce a generalization of the CTES algorithm. In particu¬ 
lar, we generalize the scaling property defining the new auxiliary variable 

^N,s = s^N = sN ^, (56) 

where s is equal to the product of one or more generic prime numbers (s = Tl^Pk)- 
We can now rescale the CTES function in Eq. Q as a function of the continuous 
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Fig. 6 Factorization of any integer in the range Nmi„ = I <N < N,„a.x = 64 by using the method 
(2) aimed at checking all the trial factors in the r ang e We exploit n — 3 interferograms 

defined by Eq. jl^ , with i — 1,2,3, where the range of values 
satisfies the condition nil and the values jcp^ are given by Eq. psl. We find that the factor 
3 of 15 corresponds to a maximum of the interferogram while the factor 5 is 

associated with a maxima of the interferogram On the other hand the factor 9 of 

63 is associated with the maximum of the interferogram 


variable 


1 M 

m=l 

Let us describe how such a rescaled function allows us to find the factors of a 
generic number N. If one of the chosen prime numbers that defines the value of s 
is a factor of N, we have solved the problem. If not, the factors of N are the integer 
values of ,, different from the prime numbers that correspond to dominant 
maxima in the rescaled CTES function given by Eq. ([57)i. 


2Ki{m — 1)-' 


Ns 

^N,s 


(57) 
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In a general analogue implementation, the auxiliary variable in Eq- < 56 1 is 
obtained by rescaling the values of the observable according to 


. _sN 

hN.s = -Of 

X ^ 


(58) 


Thereby, the range 


(max)- 

'N.S 


of trial factors covered by the variable ^n,s can 

be equivalent to the interval of trial factors covered by the variable 

by exploiting a range for the variable Of of the physical observable (9f of length 
Ao = Omax — Omin reduced by a factor s. However, at the same time, each of the values 
x=Xi, with / = 0,l,...,n—1, increases by the same factor s according to the iteration 
formula (311 where xq in Eq. (491 now reads 


Xq — ^XfnaxOmin/ (59) 

with ^0 = 1 j y/Nmin depending on the range ( |37| or ( |3^ , respectively, of trial factors 
we are considering. This implies that the maximum factorable integer Nmax is upper 
limited by xoj{somin), where Xmax is the maximum value physically allowed for the 
parameter x. 


5 Remarks 

We have described a novel analogue algorithm based on the experimental measure¬ 
ment of CTES interferograms (of ;x) depending on the continuous argument Of 
and the discrete parameter x associated with two suitable physical observables Of 
and Ox, respectively. The domain of factorable integers is determined by the range 
Omin <o^< Omax of the Continuous variable Of and by the maximum value Xmax al¬ 
lowed by the parameter x . 

The largest integeiWmm: factorable with a single CTES interferogram (of ;x) 
associated with a given value of x and defined in a certain range o,nin < Of < Omax is 
upper limited by the value {omax/ominY- For larger integers N the available physical 
spectrum of observable values Of may not be enough to cover all the trail factors ei¬ 
ther in the range 1 < < y/N or y/N< < N. Eor this reason, we have introduced 

an algorithm based on the measurement of n different interferograms 
with i = 0,1, ...,n — 1, associated with the respective values .r = Xi. In this case the 
largest number factorable is upper limited by the value Xmax/omin- We have demon¬ 
strated that the number n of necessary experimental runs scales logarithmically with 
respect to the root of the number N we want to factor. Very interestingly our method 
allows a parallel factorization of an exponential number of large integers with re¬ 
spect to the number of binary digits timax associated with Nmax- These results are very 
important in view of the optical implementation of the algorithm described in Ref. 
0 . 

Moreover, we have introduced a generalized CTES procedure defined by the scal¬ 
ing property ^N, s = sN^, where s = Fli^pk with generic prime numbers which are 
non factors of N. We have shown that it is possible to reduce for a factor s the length 
Ao = Omax —Omin of the range of values in which the function I^^’^\o^-,x) is 
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recorded. In such a case, the maximum factorable integer Nmax is upper limited by 

In Ref. Q, we describe in detail how an optical computer enables a physical com¬ 
putation of the CTES algorithm for several orders j. In such a case the values and 
X defining the analogue interferogram correspond, respectively, to the 

wavelengths A of a polychromatic source and to the unit of displacement defining the 
optical paths in a generalized Michelson interferometer IJl. Indeed, an experimental 
proof of the principle of the CTES algorithm has been performed in the case of j = 2 
and j = 3, leading to the factorization of seven-digit numbers m\. 


6 Towards a polynomial scaling in the number of physical resources 


We point out that any “classical” implementation of the algorithm described so far 
would lack exponential speed-up. Indeed, the largest number factorable Nmax is upper 
limited either by the value {omaxiominf' xojomim depending on the use of a single 
interferogram or a sequence of interferograms. 


Moreover, the accuracy in the variable ^ in Eq. (Ill 


, f; ^*3 A I ^ ^ ^max , I ^ yi 

Aq = -^Ax+ -Aof < — :^Ax+ -Aof 

x^ X ^ x^ X ^ 


depends on the given experimental indeterminations Ao^ and Ax associated with the 
measurement of the observables and Ox, respectively. This shows from Eqs. (|^ 
and (|^ that the unit x defining the CTES interferograms in Eq. ( [T^ needs to grow 
exponentially with respect to the number of bits associated with the largest number 
Nmax in the interval of integers to be factored. 

The ability to resolve the maxima associated with factors from non factors for 
larger and larger integers N Q] can be, in general, improved by implementing ’’factor¬ 
ing” interferograms of either a larger order J or with a larger number M of interfering 
terms as can be inferred by the relative curlicue functions in Eig. 1. 

Ultimately, only a “quantum” system able to encode such physical observables in 
a qubit representation in line with Shor’s algorithm would make it possible to avoid 
the requirement of exponentially large values for the observable Ox- In particular, a 
crucial role in our algorithm is played by the hyperbolic function f(^) in Eq. ([^ 
emerging from the ratio of the values Ox and and characterizing the “factoring” 
CTES interferograms in Eq. ( [T2| ). A “qubit” parallel representation of functions anal¬ 
ogous to the curlicue function may lead to novel factoring algorithms based on a 
polynomial number of resources. Moreover, multi-qubit quantum interference may 
serve as an efficient tool to distinguish factors from non factors. 
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